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The most important character of the exotic superconductor PrOs4Sbi2 is the existence of low- 
lying excitations (excitons) with a finite energy gap and it appears as the magnetic field-induced 
order above 4.5 T. We focus on the a u conduction band, which hybridizes with a Pr 4/ 2 state 
strongly, coupled to the excitons. It results in an attractive interaction between the conduction 
electrons. The symmetry of the superconducting order parameter is determined by dispersion 
relation of the exciton. For the bcc system PrOs4Sbi2, a d-wave state (k x k y +uk y k z +ui 2 k z k x , 
stabilized with broken time reversal symmetry. 
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PrOs4St>i2 is a recently found superconductor in the 
Praseodymium-based heavy- fermion system. 1 Specific 
heat measurement revealed multiple superconducting 
transition temperatures at T c i = 1.85 K and T C 2 = 1.75 
K. 2 The experiment on thermal conductivity also re- 
ported the multiple phases depending on temperature 
and external magnetic field. 3 The nuclear quadrupole 
resonance (NQR) experiment showed that there is no 
Hebel-Slichter peak in T^ 1 at the superconducting tran- 
sition temperature. 4 Very recently, zero-field /xSR mea- 
surement revealed that the superconducting state is as- 
sociated with a spontaneous magnetic field, indicating 
that the superconducting state breaks the time rever- 
sal symmetry. 5 Thus far, only phenomenological theo- 
ries have tried to account for these experimental results 
of unconventional superconductivity. 6-10 Since the mul- 
tiple superconducting phase is still under investigation, 
we do not discuss it in this letter. The main purpose of 
this letter is to present a microscopic theory considering 
electronic states specific to PrOs 4 Sbi2, and to give a sce- 
nario for the superconductivity with broken time reversal 
symmetry. 

The Pr + ion has 4/ 2 configuration in a Th point group 
crystal field. 11 It is reported that the Fermi surface of 
PrOs4Sbi2 is similar to that of LaOs 4 Sbi2 which is a ref- 
erence compound, 12 indicating well-localized 4/ 2 nature. 
This localized nature of the f-electrons is characteristic 
to PrOs 4 Sbi2 , compared with other heavy fermion super- 
conductors with itinerant /-electron nature such as in U- 
based compounds. Therefore, we consider the conduction 
electron system well-separated from the /-electrons. 

Another characteristic point of PrOs4Sbi2 is the mag- 
netic field-induced ordered phase above 4.5 T, observed 
by specific heat, 13 electric conductivity, magnetization 
and thermal expansion measurements. 14 In this ordered 
phase, high-field neutron scattering measurements re- 
vealed a small staggered moment perpendicular to the 
field. 15 Shiina and Aoki proposed that the field-induced 
order is mainly driven by a quadrupole-quadrupole inter- 
action. 16 They assumed the Ti singlet ground state and 
the r$ triplet first excited state in Pr 4/ 2 configuration 



(specifically, T 4 should be used for Th representation). 
Since an external magnetic field lifts the degenerate T$ 
triplet and stabilizes one of them, this Ti-Ps level scheme 
explains the field-induced order. In this letter, we assume 
this ri-r5 scheme. 

PrRu4Sbi2 is a reference superconductor (T c = 1.04 
K 17 ) with a Hebel-Slichter peak in Tf 1 of NQR, 18 in- 
dicating an s-wave pairing state. It is reported that 
PrRu4Sbi2 has also triplet excitations as in PrOs4Sbi2. 19 
However, a field-induced order has not been reported 
thus far in PrRu4Sbi2. This means that the crystal- field 
excitation gap to the triplet state is much smaller in 
PrOs 4 Sbi2 than in PrRu4Sbi2. The low-energy excita- 
tions play important roles in the field-induced order and 
the heavy electron mass for PrOs 4 Sbi2. The low- lying 
excitation (exciton) is expected to be the most impor- 
tant origin of the exotic superconductivity in PrOs4Sbi2, 
while this is not the case for PrRu4Sbi2. In this letter, we 
present a microscopic theory for time reversal breaking 
superconductivity mediated by excitons, specific to the 
bcc system PrOs4Sbi2. 

The wave functions for the Pr 4/ 2 state are 
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Here, Ti, and T 5 are Oh representations. The total 
angular momentum is fixed to 4, and the wave func- 
tion \J Z ) (J z = —4^4) represents the z component 
of the angular momentum. We take the basis functions 
\(j) n ) (n = 0, 1, 2, 3) for the T h system as 16 



der in the boson operators as 



Hf = A CF °in°* 



(7) 



I0o) = |ri>, 

I0n>o) = Vi-^ir^ + dirj), 
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where — 1 / y/2 < d < 1/ -\/2- As an effective Hamiltonian 
for the 4/ 2 states, we take an intersite interaction into 
account as well. 

Hf = Hqf + Hi 

3 

H CF = A CF ^ \<j> n ) (4> n | (3) 
n=l 

H *= E EE D PP X ip X iP 

a=l,2,3 <ij) /3 

Here, i?cF is a local term with a crystal field excita- 
tion gap Acf- We choose isotropic multipolc-multipolc 
interactions for H\ in order to simplify the following dis- 
cussion. In Hi, Dp/] is a coupling constant and Xp is one 
of the three-dimensional tensor operators (denoted by a) 
for a dipole, a quadrupole and an octupole, among others 
(denoted by (3). T^Uj) denotes the summation over the 
nearest-neighbor Pr sites in the bec lattice. 

First, we introduce bosonic operators such as 20 ' 21 



A E a ln a jn + >>'{a\l a )z + 43 fl ]l - 4^2) + h - C - 
(ij) L n=l 

Here, A' = A = x^Dqq holds for the only quadrupole- 
quadrupole interaction in H\. For a dipole-dipole inter- 
action, A' = —A is obtained. In general, other multipolcs 
are involved in the intersite interaction as well as the 
quadrupole, so that A' differs from A. In our theory, the 
types of intersite interaction do not alter our result. They 
just modify the constants A and A'. 

Since a\ couples with 03 by the pair creation and an- 
nihilation terms, we transform the operators as a\ 



-(a x - ia y )/\/2, a 2 



03 



+ ia y )/V2. Hf is 



then given by the decoupled bosons a a (a = x, y, z). 
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a' ko was introduced with N 



n,n'—0 



as the number of Pr sites. We can now diagonalize the 
Hamiltonian using a Bogoliubov transformation. 



with x nn'f3 = (<Pn\X^\(j) n ') and a\a n , = \<p n )(<p n >\. The Ha = E k^ ka b ka - -E* 6 k-E k ) 



bosons are subjected to X)«=o a li a n = 1 a t each site. Due 
to the finite excitation gap, bosons a n >o for excitations 
are dilute at low temperatures, and we eliminate ao for 

the ground state using a = aj = \J 1 — Y^ n =i a «a„. The 
multipole operator is rewritten up to the quadratic order 
in dilute boson operators as 



a ka = u k b ka + v fc &i 



ka 
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One of the relevant intersite interactions in PrOs4Sbi2 

is of the T 5 quadrupolar type (Xq — O xy , Xq — O yz , 

C 3 
Q 



Xq — O zx ). W 



Ofr = —{-h J n + J n J d 



(6) 



Using this quadrupole operator, we demonstrate a 
derivation of exciton dispersion. Nonzero matrix ele- 
ments for the quadrupole operators x" n ,p (/3 = Q) 
are given as x\ 0Q = -\x Q , x 2 WQ = ix Q /V2, xj QQ = 
-1XQ/V2, xf 0Q = xq/y/2, xl 0Q = xq/V2, and a; Q = 
^35(l-rf 2 ). H f is now derived up to the quadratic or- 



Ufc = V^ ( lt + 1) ' «* = ^(A*)^(^-1) 

Here, b k describes low-energy bosonic excitations. The 
dispersion relation of the exciton is given by E k = 

y/ e 2 fc — A 2 fc . There are threefold degenerate excitations 

(a = x,y,z). We note that the Hamiltonian (8) is the 
same as that for interacting spin dimer systems in which 
the field-induced order takes place. 22 

Next, we derive an effective interaction between con- 
duction electrons via exciton creation and annihilation 
processes. The conduction electron system has charac- 
teristics of both the & u and t u molecular orbitals of the 
Sbi2 cage structure. Since the & u component strongly 
couples with the 4/ 1 wave function of the 4/ 2 T 5 triplet 
at the r point, 23 we restrict ourselves to the & u con- 
duction band to discuss superconductivity. We study the 
following Hamiltonian H. 



H = H c + Hf + H c f 



(10) 



He — E 6 fc ct fc<x C kai H cf — -J c f y^ Sj-Sj 
ka i 

Here, H c is for the conduction electron system. Hf is 
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given in eq. (8). H c f represents an effective exchange 
interaction (J c f > 0) between the conduction electrons 
and the 4f 2 T 5 triplet. 24 The Ti and T 5 states do not 
couple with each other via only spin exchange. H c f is 
isotropic in 7\ symmetry. Sj denotes a spin (S = 1/2) 
operator for a conduction electron system at the i-th site. 



(11) 



Here, Ci a is an annihilation operator for the a, u electron 
at the z-th site. Si is a pscudospin (5 = 1) operator for 
the \4> n >o) triplet, which is given at each Pr site by 

S = -i(<4,a],,<4) x (a x ,a y ,a z ). (12) 

When conduction electrons excite the triplet excitons via 
H c {, they induce polarization of the pscudospin S. The 
other conduction electrons approach the polarized site 
in the next process, giving rise to an effective interaction 
between the conduction electrons. We treat H c f as a per- 
turbation and derive an effective interaction Hamiltonian 
up to the second order. 

H ' = Hc fp — w — w H * (13) 

Since the characteristic energy is Aqf ~ 8 K, 15 ' 25 the 
excitons are dilute at T < T c — 1.85 K. We consider only 
the exciton pair creation and annihilation process which 
is dominant at low temperatures. In this case, we obtain 
the following H 1 using the Bogoliubov transformation 
(9). 
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Here, s k is the Fourier-transformed spin operator of the 
conduction electrons. We neglected the energies of the 
conduction electrons in the denominator. The excitation 
has a small dispersion, 26 namely |A k \ <C Acf, and V k 
is simplified as 



Vi 
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In real space, H' can be written as 

(ij) 1 

The first term is an antiferromagnetic interaction be- 
tween the nearest- neighbor sites [along (111) directions], 
while the second term is a ferromagnetic interaction on 
the same sites. The former leads to a singlet pairing. The 
latter is a ferromagnetic short-range (repulsive s-wave) 
interaction, and does not contribute to triplet pairings. 
We study then the effective Hamiltonian for the conduc- 



tion electron system. 

H eS = 6 kc f ka c ka + ^Vj, s k - s_ k (17) 
ka k 

A simple mean-field analysis leads to the following gap 
equation for singlet pairings at low temperatures. 

A(fc) = iyV;» v , A(fc,) (is) 



Here, k , represents the even component of V k _ k > 
under k — > k or k! — » — fe' transformations. It is ex- 
pressed as 



V k-k 



V {[f s (k)f s (k')-l] + f dl (k)f dl (k') 
+ fd2(k)f d 2(k') + .f d3 (k)f d3 (k')}, (19) 
/ s (fe) = cos(fc x /2) cos(k y /2) cos(fc z /2), 
fdi(k) = sin(fc x /2)sin(fcj,/2)cos(fc z /2), 
f d 2(k) = cos(fc :E /2)sin(fc y /2)sin(fc z /2), 
fda(k) = sm(k x /2)cos(k y /2)sm(k z /2). 

For triplet pairings, we note that effective interactions 
are all repulsive. For singlet pairings, there are s-wave 
[fs(k)f s ( fe') - 1] and d-wavc [f dn (k)f dn (k')] (n = 1, 2, 3) 
channels. The d-wave channel is attractive, while the s- 
wave channel is repulsive due to the —1 term. 

Now we determine what type of superconducting state 
is favorable for the d-wave. For this purpose, we derive 
a Ginzburg-Landau free energy from the gap equation. 
The gap equation for d±, corresponding to d xy -wave, is 
written using the Matsubara frequency. 

A(fe') 



III OC > / dtt ;•■/.;; (fc')- . . . 

£7 V^ + |A(fe')l 2 

A(fe) = m.fdi(k) + m.fd2(k) + V3.fd 3 (k) 



(20) 



Here, 771, r] 2 and 773 are complex numbers representing 
order parameters for the d xy , d yz and d zx -wave, respec- 
tively. J dfl k ' means an integral over the Fermi surface. 
We assume |A(fc')| is small near T c , and expand the de- 
nominator. The third-order terms are given by 



SF 4 



dn k ,f dl (k')\A(k')\ 2 A(k' 



(21) 



= A\ m \ 2 m + 2B{\ m \ 2 + \ m \ 2 ) Vl + B( v 2 + nDni 
A = jdn v&ik'), B = jdti k >.f 2 dl {k')f 2 d2 (k') 

Here, F4 is the fourth-order term in the free energy writ- 
ten as 

F 4 oc^A(\ m \ i + \ V2 \ i + \ V3 \ i ) 

+ 2B(\ Vl \ 2 \r l2 \ 2 + |77 2 | 2 |%| 2 + M 2 M 2 ) (22) 
+ \B{n?nl + vfvl + vfv! + «.). 

The generic form of F4 for the three-dimensional rcprc- 
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sentation is given by 27 

f 4 - /3i(M 2 + \ m \ 2 + \m\ 2 ) 2 + falvl + + %T 

+ /3 3 (|r/i| 4 + h 2 | 4 + |773| 4 ). (23) 
There are the following relations between the coefficients. 



0! cx B > 0, 



02 cx Ib > 0, 



03^\(A 



35) (24) 



For a simple spherical Fermi surface, 0s > holds. In the 
most stable state, the time reversal symmetry is broken 
with the following order parameter: 27 ' 28 

A(fe) = V [fdi(k) + u.f d2 {k) + oj 2 f d3 (k)}. (25) 

Here, u> = e ±l27r//3 . This state has eight point nodes. 
There is eightfold degeneracy due to the four directions 
of the threefold axes and to the time reversal degeneracy. 

In this letter, we studied the exciton mediated su- 
perconductivity to discuss the time reversal symmetry 
breaking state realized in PrOs4Sbi2. Our idea is based 
on the fact that there exist low-lying excitations above 
the nonmagnetic Ti singlet ground state. This is a unique 
feature of PrOs4Sbi 2 among known skutterudite super- 
conductors. Our theory reveals that the superconduct- 
ing pairing symmetry is determined by the k depen- 
dence of the exciton dispersion. For a bec lattice, a three- 
component d-wave state appears with broken time rever- 
sal symmetry, which agrees with the result of the zero- 
field fiSR experiment. We point out that the most stable 
state (25) breaks the time reversal symmetry around one 
of the threefold axes due to the relative phase e ±l2 *'/ 3 be- 
tween the three components. 

A magnetic exciton mechanism was suggested as a ori- 
gin of a (i-wave superconductivity in UP d2Ai3. 29-31 In 
this case, the superconductivity coexists with antifcrro- 
magnetic order. In our model for PrOs4Sbi2, the ground 
state is a nonmagnetic Ti singlet which does not suffer 
superconductivity, and we consider excitons with a finite 
energy gap. These points are different from the model for 
UPd 2 Al 3 . 

Our theory can predict that the gapped-exciton- 
mediated superconductivity is suppressed with the in- 
crease in crystal field excitation gap Acf [see eq. (15)], 
which is the characteristic energy for the excitons. 
This can be realized by substituting Ru for Os, since 
PrRu4Sbi2 has a much higher crystal field excitation en- 
ergy gap. 19 Very recently, it has been reported that T c 
decreases from 1.85 K (T c of PrOs 4 Sbi 2 ) with the substi- 
tution of Ru for Os, and then it increases towards 1.04 K 
(T c of PrRu 4 Sbi 2 ), which is understood as a competition 
between the d-wave and s-wave. 19 The first decrease in 
T c agrees with our theory. Another way to increase the 
crystal field excitation gap Acf is by applying pressure. 
In our model, the T c of PrOs4Sbi2 should be decreased 
with pressure. 

Throughout the letter, the conduction band is re- 
stricted to the o„ band hybridizing strongly with an /- 



electron wave function of the Pr 4/ 2 state. In fact, the 
t u band electron also couples with the 4/ 2 state and ad- 
mixes with the a u on the Pr sites, which leads to an or- 
bital exchange interaction. This interaction couples the 
Pr Ti singlet ground state to the T 5 triplet, and the 
orbital exchange (or interband exchange) process is in- 
volved in the effective interaction between conduction 
electrons. The possibility of multiband superconductiv- 
ity will be studied in the future. 

The idea of gapped-exciton-mediated superconductiv- 
ity can be adopted to spin dimer systems, where a non- 
magnetic singlet ground state is realized accompanied by 
dispersive triplet excitations with a finite energy gap. It 
is essentially the same as the Pr 4/ 2 system discussed 
in this letter. Introducing a conduction electron system 
coupled to the spin dimers, we can expect superconduc- 
tivity mediated by the triplet excitons. This is one of the 
directions to search for a new superconductor. 

We are indebted to Y. Aoki, K. Kuwahara, H. Shiba 
and R. Shiina for valuable discussions. 
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